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Abstract 

• We adopt Leaver's ^ method to determine quasi normal frequencies of the Schwarzschild 

. black hole in higher (D > 10) dimensions. In D-dimensional Schwarzschild metric, when 

D increases, more and more singularities, spaced uniformly on the unit circle \r\ = 1, ap- 
proach the horizon at r = r/j = 1. Thus, a solution satisfying the outgoing wave boundary 



o 
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^ ■ condition at the horizon must be continued to some mid point and only then the continued 

fraction condition can be applied. This prescription is general and applies to all cases for 
*0 ' which, due to regular singularities on the way from the point of interest to the irregular 

singularity, Leaver's method in its original setting breaks down. We illustrate the method 
calculating gravitational vector and tensor quasinormal frequencies of the Schwarzschild 
^ ' black hole in Z) = 11 and D = 10 dimensions. 

We also give the details for the D = 9 case, considered in [2]. 
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§ ■ 1 Introduction and setup. 

o ■ 

' Our motivation to study quasinormal modes of Schwarzschild black holes in higher dimen- 

Q^. sions comes mainly from the possibility of studying the dynamics of gravitational collapse 

^ I in vacuum initiated with the work of Bizoh, Chmaj and Schmidt At the expense of go- 

^P' ing to higher {D > 5) odd dimensions, they showed how to evade Birkhoff 's theorem and 

study gravitational collapse in vacuum at radial symmetry. It was shown [21 12] that for 
^ ! D = 5, 9, the D-dimensional Schwarzschild black hole is the attractor for large initial data 

^ I and at some intermediate times the solution settling down to the Schwarzschild black hole, 

obtained from nonlinear numerical evolution, is well approximated outside the horizon by 
the least damped quasinormal mode. Therefore the precise values of fundamental quasi- 
normal frequencies of Schwarzschild black holes in odd dimensions are urgently needed, as 
they help to check the validity of the numerical code used in evolution. The reliable values 
of quasinormal frequencies of Schwarzschild black hole are available for D = 5 case |11 13] , 
but for D > 5 only results from WKB methods are published [HI 13 |H1 and it is known that 
they may be not accurate for small values of angular momentum and/or higher overtones. 
Therefore getting these values with Leaver's method of continued fractions Jj, giving the 
most precise values of quasinormal frequencies in = 4 and D = 5 dimensions, is worth- 
while. We describe below how to modify Leaver's method to obtain gravitational vector 
and tensor quasinormal frequencies of Schwarzschild black hole in L> > 10 dimensions. 
The line element of the Schwarzschild solution in D-dimensions has the form 

ds^ = A{r)dt^ - A-^{r)dr'^ - dft^^^, (1) 

with 

Air) = l-(^y~\ (2) 
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type differential equation 



wliere tlie tortoise coordinate x is defined by dx/dr = A~^{r) and the parameter s 
depends on the type of the perturbation (s = for the gravitational tensor and s = 2 for 
the gravitational vector perturbation). Eq. IS)), derived independently by Gibbons and 
Hartnoll [HI and Ishibashi and Kodama 10., generalises the well known Regge- Wheeler 
equation [H| to D-dimensions. In this setting quasi normal modes are defined as solutions 
of Q, satisfying the outgoing wave boundary conditions 

V' °° exp(ibiA:x), (4) 

with Im(/c) < 0. The corresponding values of k are called quasi normal frequencies. Eq. 
© has D — 2 regular singular points (at r = and at D — 3 roots of r^~^ = 1) and 
the irregular singular point at infinity. Leaver's method of determining quasi normal 
frequencies consists in separating boundary behavior and then transforming r into p{r) 
in such a way that the singularities at r = 1 (horizon) and at r = oo become the closest 
singularities in the p plane. In D = 4 it is accomplished with the substitution 

4,{r) = (r - l)-*V2ifcgi/=r^ ^ 

where 

oo oo / 1 \ ^ 

^(p(0) = E«nP" = E«4^ ' (6) 
n=0 n=0 V / 

and the coefficients are given by the 3-term reccurence relation 

j^r^^an+i + -fl^^an + T^^^an-i = 0, (7) 

with the initial condition uq = 1, a_i = and 7n given in Then, the quantization 
condition is the convergence of the series © on the convergence radius p = 1- The two 
linearly independent solutions of the reccurence behave as 



~ exp 



(±V-8ikn^ , (8) 



thus the minimal solution makes the series © convergent at p = 1. The discrete values 
of k for which the solution given by the initial condition qq = 1 and a_i = is a minimal 
one define quasi normal frequencies. For these values the following equation, involving an 
infinite continues fraction holds 

«» 'r- 7f- ■■■■ " 

To determine the quasi normal frequencies we truncate the infinite continued fraction in 
© at some denominator and seek for the solutions of Q which are stable with respect 
to the change of depth of this truncation. 
In general for even D > 4 the substitution 

\ / n=0 

leads to a {2{D — 3) + l)-term reccurence relation, while for odd D the substitution 
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more of the D — 3 singularities, spaced uniformly on the unit circle |r| = 1, approach 
the horizon at r = 1 and no simple transformation can move them away from the circle 
centered at the horizon and the radius corresponding to r = oo in the p plane. In the 
case of eq. (jSJ, this difficulty arises first at D = 10. (The transformation p{r) has to be 
simple enough to be easily inverted into r{p), in our case it is a homography). Therefore, 
for D > 10 the solution starting from the horizon has to be continued through some 
mid points < p < 1, laying within the convergence radius of the presently used series 
representation of the solution, and Leaver's continued fraction condition can be applied 
only if it is the irregular singularity corresponding to r = oo, which limits the convergence 
radius of the presently used series representation. 



2 The D = ll case. 

As an example, to illustrate how the above prescription works, we determine quasinormal 
frequencies of the Schwarzschild black hole in D = 11 dimensions. We choose D = 11 as, 
due to our motivation given in the introduction, we are interested in odd dimensions and 
D = 11 is the smallest odd dimension in which Leaver's method in its original setting 
breaks down for D-dimensional generalization of Regge- Wheeler potential (jSJ. In eq.® 
we substitute 

/r — 1 \ ., r — 1 

— T e^'^^(pW), P{r) = • (12) 

\r + 1 J r 

The singular points of eq.© at r = 1, e^™/^, oo are transformed into /5 = 0, 1 — l/\/2 i 
i/-v/2, 1 respectively (other singular points are placed at \p\ > 1). The singularities at 
p = 1 — l/\/2 lb limit the convergence radius of the series representation of the 



solution u (p) = J2'^=o ^nP^\ to Y 2 — \/2 0.76. We choose p^ = 1/2 (which is a regular 
point) as a mid point. We have 



oo 



u{p) = Y,anP'' = Y.~^n{p- P,T, (13) 



n=0 n=0 

where 



oo oo 

Oo = X! ^-^P^' ^1 = X! "nanPo"^- (14) 

n=0 n=l 

The coefficients fulfill 2{D — 3) + 1 = 17-term reccurence relation, which reduced to 
the 3-term one via Gauss elimination ^^1) and then inserted into @ yields quasinormal 
frequencies. 

All reccurence relations are obtained analytically in Mathematica computer algebra pack- 
age. Then all other tasks (finding di/aQ from eq. (|14j) . reduction of the (2(D — 3) + l)-term 
reccurence relation for a„ to the 3-term one via Gauss elimination and finding the roots 
of the continued fraction relation 0) are done numerically, by the program in the C 
programing language. To determine the roots of the continued fraction relation Q we 
use Newton- Raphson root searching algorithm |13j . 

The first three quasinormal frequencies for vector and tensor gravitational perturbations 
of the Schwarzschild black hole, for different values of L, are given in tabled Our values 
of fundamental frequencies for tensor modes are consistent with [0] (in [H] the values of 
fundamental frequencies for scalar field perturbations were given and the scalar field per- 
turbation obeys exactly the same equation as the tensor gravitational perturbation). In 
order to compare the values obtained from our modification of Leaver's method with the 
values given in |S] we also calculate the first three quasinormal frequencies for vector and 
tensor gravitational perturbation of the Schwarzschild black hole in = 10 dimensions. 
They are listed in table |2j Comparing with [Hj we see a perfect agreement for larger 
values of L. This makes us feel cofident in our results. However for smaller L values. 
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3 


4.4533 - 


1.0331j 


3.4147 - 


2.9417i 


1.9955 


- 3.7743i 


5.1231 - 


1.0507i 


4.2669 - 


3.0765i 


2.7018 


- 4.0946i 


4 


5.2343 - 


1.0187i 


4.4049 - 


2.9628i 


2.8490 


- 4.0060i 


5.8540 - 


1.0463i 


5.1305 - 


3.0910i 


3.4995 


- 4.4424i 


5 


6.0134 - 


1.0120i 


5.3226 - 


2.9745i 


3.7955 


- 4.3603i 


6.5849 - 


1.0435i 


5.9561 - 


3.0967i 


4.4432 


- 4.7715i 


6 


6.7881 - 


1.0097j 


6.1936 - 


2.9862i 


4.8328 


- 4.6430i 


7.3160 - 


1.0415i 


6.7587 - 


3.0994i 


5.4481 


- 4.9519i 


7 


7.5579 - 


1.0096j 


7.0340 - 


2.9976i 


5.8529 


- 4.7970i 


8.0471 - 


1.0401« 


7.5459 - 


3.1007i 


6.4039 


- 5.0333i 


8 


8.3232 - 


l.OlOSi 


7.8533 - 


3.0082i 


6.8146 


- 4.8826i 


8.7783 - 


1.0390i 


8.3226 - 


3.1014i 


7.3091 


- 5.0753i 


9 


9.0846 - 


1.0120i 


8.6576 - 


3.0177i 


7.7286 


- 4.9366i 


9.5095 - 


1.0383i 


9.0914 - 


3.1019i 


8.1784 


- 5.1000i 


10 


9.8425 - 


1. 01364 


9.4505 - 


3.0262i 


8.6084 


- 4.9742i 


10.2408 - 


1.0377j 


9.8544 - 


3.1021i 


9.0222 


- 5.1159i 


11 


10.5976 - 


1.0152i 


10.2348 - 


- 3.0337i 


9.4632 


- 5.0021i 


10.9721 - 


1.0372j 


10.6128 - 


- 3.1023i 


9.8473 


- 5.1268i 



Table 1: The first three quasinormal frequencies for vector and tensor perturbation of the 
Schwarzschild black hole m D = 11 dimensions. 



n — in 

U — lU 


vector modes 


tensor modes 


L 


n = 


= 


n = 


1 


n = 


= 2 


n = 


= 


n = 


1 


n = 


= 2 


2 


3.2334 - 


0.9603i 


2.0119 - 


2.7275i 


0.9784 - 


3.5136i 


3.9209 - 


0.9621* 


3.0410 - 


2.8514* 


1.5315 - 


3.5723* 


3 


3.9946 - 


0.9337i 


3.1233 - 


2.7391i 


1.5680 - 


3.5028i 


4.6311 - 


0.9555* 


3.9309 - 


2.8507* 


2.2219 - 


3.9757* 


4 


4.7607 - 


0.9211i 


4.0833 - 


2.7243j 


2.4839 - 


3.9023* 


5.3414 - 


0.9515* 


4.7545 - 


2.8456* 


3.1902 - 


4.4737* 


5 


5.5225 - 


0.9165j 


4.9654 - 


2.7248i 


3.6165 - 


4.2895* 


6.0519 - 


0.9489* 


5.5446 - 


2.8411* 


4.3063 - 


4.6452* 


6 


6.2781 - 


0.9158i 


5.8018 - 


2.7316i 


4.7032 - 


4.4370* 


6.7627 - 


0.9472* 


6.3149 - 


2.8377* 


5.2808 - 


4.6848* 


7 


7.0278 - 


0.9168* 


6.6097 - 


2.7400i 


5.6773 - 


4.5027* 


7.4735 - 


0.9460* 


7.0722 - 


2.8350* 


6.1763 - 


4.6985* 


8 


7.7725 - 


0.9185j 


7.3984 - 


2.7484i 


6.5828 - 


4.5412* 


8.1845 - 


0.9451* 


7.8206 - 


2.8330* 


7.0262 - 


4.7042* 


9 


8.5128 - 


0.9204i 


8.1735 - 


2.75612 


7.4452 - 


4.5675* 


8.8955 - 


0.9444* 


8.5624 - 


2.8315* 


7.8469 - 


4.7067* 


10 


9.2495 - 


0.9222j 


8.9385 - 


2.7631i 


8.2785 - 


4.5869* 


9.6066 - 


0.9439* 


9.2994 - 


2.8303* 


8.6472 - 


4.7079* 


11 


9.9832 - 


0.9240* 


9.6957 - 


2.7692i 


9.0911 - 


4.6021* 


10.3178 - 


0.9435* 


10.0327 - 


2.8293* 


9.4327 - 


4.7084* 


12 


10.7144 - 


' 0.9256i 


10.4469 - 


2.7745i 


9.8882 - 


4.6143* 


11.0290 - 


0.9432* 


10.7629 - 


2.8285* 


10.2070 - 


- 4.7085* 


13 


11.4434 - 


- 0.9270i 


11.1932 - 


2.7792i 


10.6733 - 


- 4.6243* 


11.7402 - 


0.9429* 


11.4908 - 


2.8279* 


10.9725 - 


- 4.7085* 


14 


12.1707 - 


' 0.9282i 


11.9354 - 


2.7833i 


11.4490 - 


- 4.6327* 


12.4514 - 


0.9427* 


12.2166 - 


2.8273* 


11.7310 - 


- 4.7084* 


15 


12.8963 - 


' 0.9293i 


12.6743 - 


2.7869i 


12.2170 - 


- 4.6399* 


13.1627 - 


0.9425* 


12.9409 - 


2.8269* 


12.4839 - 


- 4.7082* 



Table 2: The first three quasinormal frequencies for vector and tensor perturbation of the 
Schwarzschild black hole in L) = 10 dimensions. 

and especially for the overtones there are differences exceeding 10%. As Leaver's method 
works well both for smaller and larger values of L, we believe that the error is on the 
WKB method side (see the comments in |H]). 

3 The D = 9 case. 

Here we give the details for the D = 9 case, skipped in |2j. The substitution (fTT|) leads to 
2(D — 3) = 12-term reccurence relation 

7li)a„+i + Ti^)^^ + ... + ^(^^)an-io = 0, (15) 

with 

7^^) = 216(1 +n)(ifc- 3(1 + n)) 

= -9(24(5 + 12n)fA; + 40/c2- 252 + 147s^-12L(L + 6) -36n(9 + 13n)) 

7<i^) = 6(12(-37 + 134n)ik + 2Uk^ - 1323 + 17645^ _ 36L(L + 6) + 36n(41 - 62n)) 

ryi^) = _9(8(-289 + 292n)ik + 452^;^ - 4773 + 4312s^ - 28L(L + 6) + 12n(515 - 253n)) 

7^ = 4(6(-2186 + 1279n)ik + IbOSk"^ - 34290 + 216095^ - 36L(L + 6) + 18n(1870 - 547n)) 

7^6) ^ -(24(-3219 + 132in)ik + 6052^2 _ 265050 + 130095^^ - 36L(L + 6) + 36n(5579 - 1161n)) 

'yP = 2(12(-3109 + 986n)ik + 2080A;2 _ 167877 + 69237^^ + 144n(712 - 115n)) 

7^8) = _(48(-l007 + 260n)ik + 1888fc2 - 289215 + 1053995^ + 36n(4097 - 541n)) 

7^9) = 8(3(-849 + 185n)ik + 64^;^ - 21168 + 70565^ + 9n(1029 - 115n)) 

7^ = _2(48(-53 + Wn)ik + 32k'^ - 32463 + 101435^ + 18n(691 - 67n)) 

7^11) = 6(16(-6 + n)ik - 2463 + 7355^ + 24n(35 - 3n)) 

7(^12) ^ 9(13-2n + 7s)(13-2n-7s) 
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3 


3.5389 - 


0.8278i 


2.8438 - 


2.4808j 


1.1129 - 


3.2806i 


4.1342 - 


0.8541i 


3.5853 - 


2.5825* 


1.5262 - 


4.2850* 


4 


4.2863 - 


0.8180i 


3.7573 - 


2.4494i 


2.3546 - 


3.9155j 


4.8200 - 


0.8506* 


4.3624 - 


2.5661* 


3.2201 - 


4.3333* 


5 


5.0260 - 


0.8161i 


4.5953 - 


2.4469i 


3.5813 - 


4.0450i 


5.5061 - 


0.8483* 


5.1123 - 


2.5559* 


4.2023 - 


4.3007* 


6 


5.7578 - 


0.8171i 


5.3913 - 


2.4524j 


4.5724 - 


4.0828i 


6.1926 - 


0.8469* 


5.8463 - 


2.5491* 


5.0773 - 


4.2796* 


7 


6.4826 - 


0.8191i 


6.1617- 


2.4597j 


5.4646 - 


4.1045i 


6.8792 - 


0.8458* 


6.5699 - 


2.5443* 


5.8994 - 


4.2651* 


8 


7.2018 - 


0.82142 


6.9152 - 


2.4669i 


6.3037 - 


4.1198i 


7.5660 - 


0.8451* 


7.2863 - 


2.5409* 


6.6898 - 


4.2547* 


9 


7.9165 - 


0.8236i 


7.6569 - 


2.4735i 


7.1097 - 


4.1317i 


8.2529 - 


0.8445* 


7.9975 - 


2.5383* 


7.4592 - 


4.2470* 


10 


8.6275 - 


0.8255i 


8.3899 - 


2.4793j 


7.8933 - 


4.1413i 


8.9399 - 


0.8441* 


8.7048 - 


2.5363* 


8.2136 - 


4.2411* 


11 


9.3354 - 


0.8273i 


9.1161 - 


2.4844j 


8.6607 - 


4.1492i 


9.6269 - 


0.8438* 


9.4091 - 


2.5348* 


8.9571 - 


4.2365* 


12 


10.0409 - 


- 0.8288j 


9.8370 - 


2.4887i 


9.4161 - 


4.1558i 


10.3140 - 


- 0.8435* 


10.1111 - 


- 2.5336* 


9.6922 - 


4.2328* 


13 


10.7442 - 


- 0.8301j 


10.5537 - 


- 2.4925i 


10.1620 - 


- 4.1614i 


11.0011 - 


' 0.8433* 


10.8112 - 


- 2.5325* 


10.4207 - 


- 4.2298* 


14 


11.4458 - 


- 0.8313j 


11.2670 - 


- 2.4957i 


10.9004 - 


- 4.1662i 


11.6883 - 


' 0.8431* 


11.5098 - 


- 2.5317* 


11.1439 - 


- 4.2274* 


15 


12.1460 - 


- 0.8323j 


11.9774 - 


- 2.4986i 


11.6328 - 


- 4.1703i 


12.3755 - 


- 0.8430* 


12.2070 - 


- 2.5310* 


11.8628 - 


- 4.2254* 



Table 3: The first three quasinormal frequencies for vector and tensor gravitational perturbation of 
the Schwarzschild black hole in D = 9 dimensions. 

Gauss elimination and insertion into Q yields again quasinormal frequencies. The 
first three quasinormal frequencies for vector and tensor gravitational perturbation of the 
Schwarzschild black hole, for different values of L, are given in table El 

4 Summary. 

We have shown how to use Leaver's ^ method of continued fraction in the case of a 
number of regular singular points, which set lower bounds on the convergence radius of 
the series representation of a solution than irregular singular point. This prescription is 
general and together with Leaver's method makes a powerful tool in determination of 
quasinormal frequencies for wave equations and resonances in quantum mechanics. 
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